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O ! 

We discuss how an extended foreground of the cosmic microwave background (CMB) can account 
£\) ■ for the anomalies in the low multipoles of the CMB anisotropies. The distortion needed to account for 

the anomalies is consistent with a cold spot with the spatial geometry of the Local Supercluster (LSC) 



and a temperature quadrupole of order ATf ~ 50 fiK 2 . If this hypothetic foreground is subtracted 
from the CMB data, the amplitude of the quadrupole {£ — 2) is substantially increased, and the 
statistically improbable alignment of the quadrupole with the octopole (I — 3) is substantially 
weakened, increasing dramatically the likelihood of the "cleaned" maps. By placing the foreground 
on random locations and then computing the likelihood of the cleaned maps we can estimate the 
most likely place for this foreground. Although the 1-year WMAP data clearly points the location 
of this hypothetical foreground to the LSC or its specular image (i.e., the vicinity of the poles of the 
cosmic dipole axis), the three-year data seems to point to these locations as well as the north ecliptic 
1 pole. We show that this is consistent with the symmetries of the cosmic quadrupole. We also discuss 

' a possible mechanism that could have generated this foreground: the thermal Sunyaev-Zeldovich 

, effect caused by hot electrons in the LSC. We argue that the temperature and density of the hot 

I/-} ■ gas which are necessary to generate such an effect, though in the upper end of the expected range 

of values, are marginally consistent with present observations of the X-ray background of spectral 
\Q ' distortions of the CMB. 
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PACS numbers: 98.80.-k, 98.65.Dx, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 



The cosmic microwave background (CMB) anisotropies have been measured with very high accuracy by WMAP 
0,0. Such a barrage of new data seldom brings only confirmation of known theories and mechanisms, and WMAP is 
no exception: lack of higher correlations 0] and some curious correlations between large-scale anisotropies are some 
of the most intriguing questions that have been raised by the WMAP data. In particular, two problems have been 
pointed with present CMB observations, using a wide variety of data, methods, maps and sky-cuts: first, that the 
quadrupole (£ = 2) has a lower-than-expected signal H, M, M. Second, that the quadrupole and octopole (£ = 3) 
present an unexpectedly high degree of alignment 0, 171 IsL I lH ITU Il2l IT^ | . The combined statistics of these effects 
implies that our CMB sky is only within the ~ 0,01% of randomly-generated maps with such anomalous quadrupole 
and octopole. 

It is important to note, first, that these large-angle anomalies were already present in the COBE data ^4|, and they 
were confirmed by WMAP both in the 1-year and in the 3-year data releases 0, Q • Second, even though the amplitude 
of the octopole C3 has increased in the three- year compared to the first-year WMAP data, the statistical relevance of 
the deviant multipoles £ = 2 and £ = 3 has remained practically unchanged in the newly released three-year WMAP 
data , whereas the outliers of the first-year data around £ w 20, £ ~ 40 and £ w 200 have either disappeared or 
become much less of a source of concern in the 3-year data. 
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These large-angle anomalies have motivated many ingenious explanations, such as compact topologies 0, 0, 0] , a 
broken or supressed spectrum at large scales 01 an d oscillations superimposed on the primordial spectrum of density 
fluctuations When the low value of the quadrupole is combined with the alignment of the directions defined 

by these two multipoles (their "normal planes" - see pil 171 l&L M. Hoi [ill l20|). the overall chance of such a statistical 
fluctuation is approximately 0.005% - 0.02% depending on the map and on the mask - i.e., only about one in 10000 
randomly generated models have a lower C'2 and a more aligned quadrupole and octopole than the observed CMB 
sky. 

As first noted by de Oliveira-Costa et al. |4|, the directions preferred by the quadrupole and the octopole point 
roughly towards the Virgo cluster — which is in the general vicinity of the dipole and the equinox, and has been 
dubbed the "axis of evil" 9J. These large-scale anisotropies appear when one compares the northern/western galactic 
hemisphere (where Virgo and most of the Local Supercluster lie) with the southern/eastern hemisphere as well. 
Unusually high correlations with the ecliptic have also been reported 

In this paper we investigate whether these large-scale anomalies can be due to extended (large-angle) foregrounds 
that have so far escaped detection. In particular, we re-examine a speculation by Abramo & Sodre [2£j that the 
explanation for the observed properties of the quadrupole and octopole is a diffuse, large-angle CMB foreground 
spatially correlated with the region of the sky occupied by the local supercluster (LSC) - which is a spot roughly 
50° x 30° centered around Virgo, at (l,b) = (284°, 74°). The peak temperature of this foreground would have to be 
of order ~ 20/iK, with a root-mean-square average temperature of order ~ 8/1K. Such a foreground would have the 
correct geometry to affect the quadrupole and octopole in a positive way, erasing their alignments and significantly 
increasing the amplitude of the quadrupole. One of the physical processes that could produce such a foreground is 
the thermal Sunyaev-Zeldovich effect (SZe) due to hot electrons in the intra-supercluster (ISC) medium j2^. For the 
range of frequencies observed by WMAP and COBE, the ISC gas causes an apparent decrease in the temperature 
of the CMB photons in the direction of the LSC. We have estimated, using a simplified model, that a temperature 
distortion as high as \ATi = 2\ rm s ~ 8 /iK is marginally consistent with present constraints on the ISC medium and 
with spectral distortions of the CMB. 

The possibility that foregrounds could be responsible for the alignments was also noted in j2|J, and two recently 
related explanations were proposed by Rakic, Rasanen and Schwarz |27| and by Inoue and Silk |28j . Rakic et al. 
studied the Rees-Sciama effect due to a nearby large-scale structure (or structures), and found that the ensuing 
foreground could be as high as ~ 30/iK, but that the phases were not right to eliminate the alignments and the low 
quadrupole. 

Inoue and Silk, on the other hand, speculate that the non-gaussian cold spot in the southern galactic hemisphere 
[2^ could be caused by voids in the nearby large-scale structure [28|. They conclude that the contribution from 
compensated pairs of voids would have the right phases to account for the low amplitude of the quadrupole and the 
for the unusual quadrupole-octopole alignment. 

Notice that, because the quadrupole is even under parity transformations n — > —ft, any given pattern has the same 
quadrupole components as its specular image. Since the effect discussed here relies mostly on a distortion of the 
quadrupole (which is both low-amplitude and has fewer phases than the octopole), this partly explains the apparent 
equivalence between the foreground proposed by Abramo & Sodre and foregrounds located in the southern galactic 
hemisphere, such as those proposed by Inoue & Silk. Furthermore, because most of the power of the observed CMB 
quadrupole lies in its m = and m = ±2 components, another probable spot for an extended foreground corresponds 
to the location of Virgo, rotated 180° around the galactic poles axis - that is, (I, b) « (100°, 70°). 

By placing one of these hypothetical foregrounds (henceforth HFg) on random locations, removing it from a CMB 
map and then computing the likelihood of the "cleaned" map, we have been able to test the randomness of the spatial 
correlation with the four dual points described above. We find that, in all maps and in all foreground models tested, 
the most probable places for them are indeed either the vicinities of ulsc, ^n^sc or near the ecliptic poles. 

We will show that, of all possible locations for this HFg, the LSC and its dual points produce the most significant 
improvements in the likelihoods of the CMB maps, by both increasing the (too low) level of the quadrupole Ci and by 
weakening the (too high) quadrupole-octopole alignment. This can be achieved with HFg's whose rms temperatures 
lie in the range 5 — 15/xK and quadrupoles C.f 9rd ~ 50 — 120/iK 2 , depending on the model. 

We have analysed the WMAP (1-year and 3-year) Internal Linear Combination maps 0,0] (henceforth ILC), the 
map of Tegmark et al. [HEl (henceforth TOH), as well as the co-added maps based on 1-year and 3-year WMAP 
data (henceforth Coadded) For the ILC and Coadded maps we use the Kp2 maks, and for the TOH map we use the 
masks M0 and M6 described in |2l|. In all cases the low value of the quadrupole and the alignments are robust, and 
removal of the HFg leads to dramatic increases in the likelihoods of the CMB maps. This strongly argues in favour 
of still unknown diffuse, large-angle structures around the dipole axis that may be affecting the CMB. 

This paper is organized as follows: In Sec. II we summarize the multipole vector formalism and the statistics of 
alignments for CMB maps. In Sec. Ill we present a model of the HFg based on the LSC, and argue that it may be 
due to the Sunyaev-Zeldovich effect caused by hot gas in the intra-supercluster medium. In Sec. IV we consider the 
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spatial location of the foreground, and show that the association with the dipole axis is not an accident. We conclude 
in Sec. V. 



II. LOW QUADRUPOLE AND ALIGNMENTS 

Katz and Weeks j2£| have described^jn a beautiful paper, how to compute all multipole vectors given the spherical 
harmonic components ae m - see also [23|. The multipole vectors, introduced to CMB data analysis by Copi et al. 
0, are essentially eigenvectors of a simple set of algebraic equation whose parameters are the multipole components. 
Very similar computations were conducted using other (usually numeric) methods in 0, Q to find these vectors. 
The idea, which goes back to J. C. Maxwell in the XIX th century, is that the multipole decomposition of a field on 
S 2 implies that for each moment £ there are £ eigenvectors of norm unity, n^ :P K The bottom line of the multipole 
vector analysis is that the expansion in spherical harmonics is equivalent to an expansion in multipole vectors: 



AIK ^ } V) a„„Y„„:0.^ />, If" 
-I p=l 



Y em (6, <p) = D t ]J ™ {e ' p) ■ H6, 0) - Z t -i(9, <f ) , (1) 



where Zi-\ just subtracts the residual £' < £ total angular momentum parts of the product expansion, and is irrelevant 
to our analysis - see j2fl| for an enhanced discussion of the multipole vector expansion. 

Notice that, contrary to the Ce's, which are always positive-definite, the Di's can be either negative or positive. This 
means that the multipole vectors n^ e ' pS) define only directions |2C|. hence they are in fact "vectors without arrowheads" . 
It can also be seen from the expansion above that, whenever using the multipole vectors to test for alignments, it is 
irrelevant what the amplitudes of the multipoles are - just their (complex) phases matter, of which there are £ for 
each multipole. 

Starting with these £ multipole vectors one can also construct £{£—!)/ 2 normal vectors - or normal planes. Therefore, 
for £ = 2 there are 2 multipole vectors (tt/ 2 ' 1 ) and n' 2 ' 2 * 1 ) and only one normal plane (u^ 2,1 ) = n^ 2 '^ x r^ 2,2 )); for £ — 3 
there are 3 multipole vectors and 3 normal planes; and so forth. Notice that, because the multipole vectors are not 
necessarily orthogonal, the normal vectors need not be (and generally are not) of norm unity. 

We can therefore check for "alignments" between either the multipole vector themselves, or between the normal 
planes. Two widely used tests that check for alignments of the quadrupole and octopole normal planes are the S 
statistic: 

S = I |^(2,l) . ^(3,1) | + |^(2,1) . 0(3,2) | + |^(2,1) . -(3,3) | ^ (2) 
o 

and the D statistic, which is analogous to S but disregards the norm of the normal vectors: 

= Ilrf,(2,l) . 0(3,1)1 + 10(2,1) . 0(3,2), x |,~(2,1) . ,-,(3,3) 
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+ |0(2,D .0(3,3^ _ (3) 



It can be easily seen that both S and D lie within the interval (0, 1). In what follows we will use mostly the S statistic, 
since the D statistic disregards the norm of the normal vector and therefore throws away some information about the 
system. For the higher multipoles this may not be much of an issue, but the quadrupole has only 2 complex phases 
(4 effective degrees of freedom) and we would like to retain as much of that phase information as possible. 

One can easily compute the probability distribution functions (P.D.F.'s) for S and D using randomly-generated 
("mock") maps, either by simulating maps, computing the harmonic components and then the vectors and alignments, 
or by directly simulating the harmonic components, which are Gaussian random numbers with dispersion given by 
(|ci£m| 2 ) = Cg. Because the alignments do not depend on the Ci's, but only on the phases between the ai m 's, one 
need not adjust the level of the Cg for each £ - all that is needed to test for the alignments in mock maps is the fact 
that the phases are random within each multipole. 

In Fig. 1 we show the normalized P.D.F.'s for the S and D tests that were computed using 300.000 mock maps. 
As mentioned above, in what follows we will use mostly the S statistic, but we note that the D statistic gives similar 
results. For a complete treatment of statistical tests of isotropy, their P.D.F.'s, and an assessment of other sources of 
error see 1241. 



Statistics of large-angle anisotropics and the low quadrupole 

Given a CMB map, the harmonic components can be easily extracted (we use HEALPix [3(j), and the multipole 
vectors and their statistics can be easily computed. There are several maps to choose from, the most well-known 
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FIG. 1: Normalized P.D.F.'s for the S statistic (left panel) and D statistic (right panel), found by simulating 3 x 10 5 mock 
maps, binned in intervals of 0.01. 
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TABLE I: Quadrupoles and alignments of CMB maps with and without the hypothetical LSC foreground subtracted. Shown 
are the Coadded and ILC maps with the 3-year KP2 mask 0,0, and the TOH map with and without the mask described in 
[iH ] (based on first year WMAP data). P_(C 2 ) is the probability that a random map has quadrupole as low as C2, and P+{S) 
is the probability that a random map has a quadrupole-octopole alignment as high as S. Also shown (last column) is the 
unbiased joint probability Prot = 16 x P_(C 2 ) x P+(C 2 ) x P-(S) x P+(S), which estimates the likelihood that a random map 
has an anomalous (too high or too low) quadrupole and an anomalous (too high or too low) quadrupole-octopole alignment. 
In all cases, removal of the foreground leads to an improvement of about two orders of magnitude in P-rot- 

being the Coadded, ILC, LILC, TOH and the Q-, V- and W-band frequency maps. For all except the TOH map 
(which is already cleaned) we have use the KP2 mask based on three- year WMAP data pj. 

It must be noted that the relativistic Doppler correction to the quadrupole is an important factor that must be 
subtracted from the maps, since it corresponds to a non-primary source of the quadrupole . 

In Table I we present the quadrupoles and their alignments with the octopoles, for a few CMB maps. The Coadded 
and ILC maps use 3-year WMAP data, while the TOH map uses the first-year WMAP data only. The alignments are 
robust in all maps, as has been noted by [8l IsL IToj]. In the next section we will construct a model foreground based 
on the LSC, and Table I presents the statistical effect of the subtraction of this LSC-shaped fore ground. Notice that 
the relative error for the probabilities can be estimated as AP(X)/P(X) ~ 1/^/300.000 x P(X). 

III. HYPOTHETICAL FOREGROUND 

As first noted by de Oliveira-Costa et al. 0, ^J, both the quadrupole and the octopole seem to be aligned on 
the plane defined by the direction (l,b) w (250°, 60°), which is quite close to the Virgo cluster. This motivated the 
proposal of Abramo & Sodre [2{| , who speculated that the low-£ anomalies of the CMB could be explained by the 
Sunyaev-Zeldovich effect caused by hot electrons in the intra-supercluster medium of the LSC. For the frequency 
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FIG. 2: Projected volume of the LSC on a Mollweide projection - the left edge corresponding to / = 0°, the middle meridian 
to I — 180° and the horizontal line to b = 0° in galactic coordinates. Virgo is at the center of the LSC, at (l,b) « (284°, 74°). 

channels observed by COBE and WMAP, this effect would cause a cold spot with the shape and location of the 
supercluster superimposed on the primary CMB data. 

We will show next that there are strong indications that an LSC-related foreground (or some structure diametrically 
opposite to the LSC) is distorting the observed CMB sky. Whether or not that foreground is caused by the SZe [25j, a 
void [23 , some other mechanism such as the Rees-Sciama (or Integrated Sachs- Wolfe) effect [2?J , or even a combination 
of those, remains to be seen. 



The morphology of the LSC is relatively well known [3j]]: it is a flattened collection of groups and clouds of galaxies 
centered at the Virgo Cluster, which contains ^20% of its bright galaxies. The Local Group is dynamically linked 
to the LSC, and lies ~15 Mpc from Virgo, at the border of the LSC. Notice that the LSC itself is not a virialized 
structure, hence the gas in its midst is not necessarily in equilibrium. 

Since we are interested in an analytic approach at this point, a radical simplification will be made, approximating 
the shape of the LSC by an oblate spheroid of maximal radius 20 Mpc with approximate axial ratios 6:2:1 [3l| . 
Therefore, our simple model assumes that the LSC is a collection of objects (clouds, groups and the Virgo cluster) 
which are distributed smoothly across the spheroid. The Sun stands at the margin of the spheroid (which looks like 
a flattened rugby ball), approximately 15 Mpc away from Virgo. 

In our foreground model we assume that the intensity of the temperature decrement is proportional to the volume 
of the LSC, projected along the line of sight (i.e., the surface density.) This must be roughly correct, whatever the 
source of the hypothetical foreground, if it is indeed correlated with a diffuse structure such as the LSC. It is trivial to 
compute the surface density, given the shape and orientation of the LSC, and the result of this projection can be seen 
in Fig. 2 for an arbitrary (but constant) density. In that figure it can also be seen that our LSC model's projection 
on our sky is a spot of roughly 50° x 30°. 

The oblate spheroid, in a conveniently rotated frame, is defined by the surface: 



where A is the major axis (along the direction z'), and b and c are the ratios of the minor axes to the major axis. 

The LSC parameters are A ps 20 Mpc, b « 3 and c sa 6. With these values the semi-major axes are B = A/b 6.7 
Mpc and C = A/c « 3.3 Mpc. Assuming that the Sun is located at a distance R under the z'-axis of the spheroid, 
the distance to the surface of the spheroid along lines-of-sight emanating from the Sun are given by: 



A. Shape and location of the Local Supercluster 



{bx'f + {cy'f + z- 
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FIG. 3: Angular power spectra of the LSC foreground model (left panel), where we use a — 1 - see Eq. Hill : and of the Disk 
(uniform-temperature) foreground model, with angular diameters of 30° (diamonds), 50° (stars) and 70° (squares.) 



Obviously, if the density is uniform inside the spheroid then the surface density will be proportional to rg. 

We can use the surface density of the LSC as the blueprint for a foreground, and therefore consider the temperature 
decrement caused by the foreground to be proportional to r$- In Fig. 3 (left panel) we show the spectrum of 
anisotropics for such a model, where the proportionality constant is set by assuming that the effect is caused by 
scattering of the CMB photons by hot electrons in the ISC medium - see below. Notice that the quadrupole is 
substantially higher than the other multipoles, because the foreground's temperature is not constant as the line of 
sight moves away from the center of the LSC. This sorts out the quadrupole as the biggest contribution to the power 
spectrum, and in fact the RMS temperature fluctuation over the whole sky is well approximated by the quadrupole. 

Given our ignorance about the existence, shape and form of this hypothetical foreground, we could equally well 
assume, following Inoue & Silk |28|. that the temperature of the foreground is approximately uniform, falling quickly 
to zero away from the center of the spot. In Fig. 3 (right panel) we show the anisotropy spectra of three such spots - 
with angular diameters of 30°, 50° and 70°. The fact that the temperature is uniform inside the spot means that the 
low multipoles contribute more evenly to the foreground. Hence, for a given RMS temperature fluctuation there is 
less power in the quadrupole in the uniform-temperature foreground relative to the varying-temperature foreground. 
Moreover, as higher multipoles are more important in the uniform-temperature foreground, if that is the case then 
it may be possible to find independent corroborating evidence by searching for anomalous alignments in the higher 
multipoles (e.g., I — 4 and I — 5) as well. This will be analysed in Sec. IV. 



B. Foreground model: Sunyaev-Zeldovich effect in the LSC 



The SZe is caused by the inverse Compton scattering of CMB photons by hot electrons in the intra-cluster medium 
|32j . It is a nonthermal, frequency-dependent effect: the upscattering causes an incident blackbody spectrum of 
photons to become distorted in such a way that the resulting higher abundance of high-energy photons is compensated 
by a shortage of low-energy photons. The spectral distortion is given by |3^ : 



AT(0, pi!/) fa x ( 



(B,tp) (zcotanh| - 4) , (6) 



where T$ — 2.726K is the temperature of the CMB, y is the comptonization parameter in the direction n{6, if) and 
x = hv/kBTo. 

The frequency at which photons are neither depleted nor overproduced is vq = 218 GHz [3j| — COBE/DMR 
and WMAP work in the range 20-90 GHz. For frequencies below Vq the effect is a nearly uniform reduction in the 
temperature of the photons, AT/T « — 2y, and for frequencies above that the effect is the opposite. This means that 
measurements over a range of frequencies around vq (such as PLANCK 's LFI and HFI can pick up the signal of 
the SZe and distinguish it from the primary anisotropies. 

The comptonization parameter y measures an optical depth for the CMB photons created by the hot electrons, and 
its value is given by the product of the Thomson cross-section ot — 6.65 x 10~ 25 cm 2 times the temperature-averaged 
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density of photons along the line of sight 33] : 



kT 

V = I <?t \n e dl , (7) 



where T e is the electron temperature, m e is the electron mass and dl = dl(0,<fr) is the line-of-sight distance element 
along the direction (8, <j)). 

The SZe has been observed over the past few years in many clusters, but its weak strength means that it could 
only be detected in the central parts of clusters, where column densities of hot gas are sufficiently high [331 l34l | . It is 
evident that some amount of SZ will take place also in the LSC, but the question is, how much? The answer depends 
on the gas density in the ISC medium, its temperature distribution, the morphology of the LSC and our position 
inside it. 

Although the morphology of the LSC as traced by galaxies is well known, the density and temperature distribution 
of the gas of the ISC medium are not. Unfortunately, X-ray and microwave observations have not yet reached the 
level of sensitivity required to detect directly the very smooth, diffuse columns of hot gas in the outer regions of 
clusters. It seems, however, obvious that there must be a great amount of ionized gas in the ISC medium, among 
other reasons because the absence of observations of the Gunn-Peterson effect indicates that most of the ISC hydrogen 
must be ionized. The gas is thought to have been shock-heated at the time of galaxy formation, and now it is probably 
distributed in many phases, including filaments and a more homogeneous component [35l l36l I37T]. Phillips, Ostriker 
and Cen [38] have constrained the amount of gas in filaments using numerical simulations and the X-ray background, 
and argued that this "warm-hot" (kT « 100 eV - 10 keV) gas can account for only 5-15 % of the "missing baryons". 
More recently, Nicastro et al. showed that this fraction could be as high as 27% j3^. It is therefore quite possible 
that much of this gas is in the ISC medium. So, the questions now are: how hot is this ionized gas, and how is it 
distributed? 

Hogan was the first to propose that superclusters (and the LSC) could impact the CMB anisotropics through the 
thermal and kinetic Sunyaev-Zeldovich effects Hjj. Molnar and Birkinshaw used HEAO 1 A2 |il| and COBE DMR 
data to analyze the Shapley supercluster and found no evidence of hot (> 10 7 K) gas in the ISC medium [42^. Boughn 
|43|. on the other hand, used the HEAO 1 A2 X-ray map and a simple "pillbox" model of nearly constant electron 
density in the LSC to argue that the SZe could be as high as \ST\ ~ (17 ± 5) /xK — although he assumed a gas 
temperature in the high end of the range 10 5 -10 8 K. Kneissl et al. [44[ did study the correlation of COBE DMR and 
ROSAT X-ray data away from the galactic plane, but it is not clear that the X-ray data has enough sensitivity to 
detect the diffuse hot gas of the LSC, and, in any case, the authors analyzed a region which misses a large chunk of 
the LSC. 

Much work has been done to study the impact of the SZe from distant clusters on the CMB (see, e.g. 
It has been found that the largest contribution to the angular power spectrum from the SZe comes from the most 
massive clusters (M ~ 10 15 /i _1 M©), at scales i ~ 3000, with amplitudes i(i + l)C e /2ir as 10 - 100 [iK 2 . 

The overall number of free electrons in the LSC can be estimated given its gas fraction and mass: 

N e = , (8) 

[i e m v 

where M^sc is the LSC mass, f g is the gas fraction, fi e is the molecular weight per electron and m p is the proton 
mass. We may assume that the mass of the LSC is ~ 7 x 10 15 Mq 49]. Assuming that the Hydrogen is fully ionized 
and that the helium mass fraction is Y — 0.24, then fi e — 1/(1 — Y/2) ~ 1.14. 

The gas fraction is not very well known, but X-ray observations of clusters indicate that f g rj O.O6/1 -3 / 2 [53. The 
fraction could be different in the ISC medium, but we will assume for simplicity that the fraction in clusters is the 
same as outside. Using h — 0.7 we get finally that the total number of electrons in the LSC should be of order 
N e ~ 7 x 10 71 . 

A very simple model for the LSC gas is to neglect the gas phases and the substructures inside the supercluster, 
and to assume that the gas is uniformly distributed inside the LSC. Hence the total volume occupied by the gas is 
Vlsc = 47r/3 x ABC and therefore the average density of electrons in the LSC is approximately: 

1.3 X 10~ 5 cirr 3 . (9) 



V L 



sc 



On the other hand, the X-ray background is also an important constraint on the density and temperature of the 
ISC medium. A compilation of observations [20] gives a background flux for energies hv ~ 2 keV of approximately 
10 -25 erg s _1 cm~ 2 sr _1 Hz -1 over the whole sky. The expected flux at this energy due to thermal bremsstrahhmg 
emission from LSC gas is ~ 5 x 10~ 26 erg s" 1 cm~ 2 sr _1 Hz -1 . Since the X-ray flux is proportional to the square 
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of the electron density, if the gas temperature is indeed 2 keV, then the upper bound for the electronic density is of 
order n e « 5 x 10 -5 cm~ 3 . This corresponds to a collapse factor of only a few. 

We can estimate the comptonization parameter assuming a constant electron density across the LSC. If the gas 
has an average temperature of 2 keV then (fcT e )/m e c 2 ~ 0.004, and with a maximum line-of-sight distance (in the 
direction of Virgo) of ~ 30Mpc we obtain that the comptonization parameter is at most: 

kT n 

Ay~v T ( \) x n e x 30Mpc ~ — —2 x lO" 5 . (10) 

jn e c 5. 10 _s cm J 

The comptonization parameter can be exactly computed from Eq. J7J for our oblate spheroid model, if we assume 
that the density and temperature of the hot gas is uniform inside the LSC, and zero outside it. In that case the 
comptonization parameter is proportional to the projected distance to the surface of the spheroid, Eq. The 
resulting angular power spectrum for the SZe of the LSC in this model is given in Fig. 3. The amplitude of the SZe 
quadrupole is: 



Aff = ^-C 2 «60aVK 2 , (11) 

_ n e (kT e ) 

a ~ 5 x 10" 5 cm- 3 X 2 keV ' 

This level of temperature distortion is consistent with the COBE FIRAS limit on deviations from the blackbody 
spectrum on large angular scales |50j . 

Recently, Dolag et al. have studied the imprint of the local superclusters on the CMB via the SZe. They have 
used constrained simulations to study these local structures, and their conclusion was that the thermal SZe is too 
small, by at least one order of magnitude, to affect either the quadrupole or the octopole. It is clear that the SZe 
model presented here cannot account for the anomalous quadrupole and octopole if the density and temperature of 
the gas in the ISC medium lie near the conventionally accepted limits, which are derived in part from observations 
and in part from simulations such as those done by Dolag et al., which include baryons, gas flows and feedback 
mechanisms. Therefore, if our model is correct, then cither the mechanisms that endow the ISC medium with gas are 
still not entirely understood, or the global parameters and initial conditions have to be changed (which is less likely.) 
Another possibility is that the SZe is not the only source of foreground, in which case a combination of foregrounds 
(all spatially correlated so the effects add up) is responsible for the distortions. In fact, such a combination of effects 
is not unlikely, since all local structures are composed of multi-phase gas and other foreground sources. 



IV. GEOMETRY AND LOCATION OF THE FOREGROUND 



In Table I we showed the effect of subtracting the HFg based on the projected volume of the LSC. This HFg 
peaks at the center of the LSC and falls steadily as the line of sight moves from the center of the LSC. When the 
HFg is subtracted from a full CMB sky, both the quadrupole is enhanced and the quadrupole-octopole alignment is 
weakened. This means that, whatever the source of the HFg in the LSC (or its specular image), the geometry of the 
LSC is such that the foreground's phases for the quadrupole and octopole work in the direction of correcting for the 
low quadrupole and for the high quadrupole-octopole alignment of the CMB maps. 

We can test the spatial location, shape and orientation of this HFg, and check whether these are indeed correlated 
with the properties of the LSC, or if the improvement in the quadrupole level and in the alignments are just flukes 
that could have happened whatever the location of the HFg. This can be done by rotating the foreground maps by 
arbitrary Euler angles, and then computing the effect of subtracting them from the CMB maps. 

We employ an unbiased measure of the likelihood of the maps defined by: 

P Tot = 2 4 P + (C 2 )P-(C 2 )P + (S)P-(S) = 2 4 P+(C 2 )[1 - P + (C 2 )]P + (S)[1 - P + (S)} . (12) 

This estimator is maximal (Prot = 1) for a map whose quadrupole C 2 and alignment S arc equal to their respective 
expectation values - in which case P-{C 2 ) = P+(C 2 ) = P-(S) = P+(S) = 1/2. 

We will consider two HFg models: the LSC model described by Eq. JSJ (hereafter HFg-LSC), and a model inspired 
by the proposal by Inoue & Silk [23 - a nearly homogeneous disk with a diameter of 50° (hereafter HFg-D.) It should 
be noticed that the HFg-D model must have a higher temperature distortion in order to cause the same order-of- 
magnitude effect in the likelihoods, compared to the HFg-LSC model. The reason is simple, and can be inferred from 
Fig. 4: since the HFg-D model has nearly homogeneous temperature, all low multipoles have similar amplitudes. 
The temperature in the HFg-LSC model, on the other hand, falls steadily from the center of the LSC, and this 
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FIG. 4: Locations of the HFg-LSC model applied to the "mask 0" TOH map [ToT.l2ll ]. The largest, darkest spot marks the top 
foreground as measured by the highest values of Prot obtained after 3000 random rotations of the putative foreground. The 
top 5%, 10%, 15% and 33% foregrounds are indicated by the progressively smaller and lighter spots. The smallest, lightest 
spots mark the bottom 67% locations. Here and in Figs. 6-16, instead of the Mollweide projection, we use a simple map 
(6,tp) -» (6,tp sin6). 




FIG. 5: Same as before, in the case of the HFG-D model applied to the "mask 0" TOH map [Til. l2ll]. 

angle dependence coincides roughly with the angular dependence of the quadrupole. Therefore, the amplitude of the 
temperature distortion in the HFg-D model must be higher than in the HFg-LSC model in order to get the same 
level of the quadrupole. So, whereas in the HFg-LSC model the peak temperature distortion is of order — 90^K, with 
a quadrupole C2 ~ 60/zK 2 , for the HFg-D model to get the same quadrupole the temperature of the 50° x 50° spot 
must be Tmsk ~ — 130/xK. 

For each model we rotate the putative foreground by random Euler angles and subtract it from a CMB map to 
obtain the "cleaned" map. We then look for the corrected maps with highest likelihoods, as measured by the largest 
value of the unbiased probability Pxot- We used 2000 random rotations for each foreground model and for each map 
of Table I (all with the KP2 mask applied.) 

The results are shown in Figs. 5-16, where each spot marks the location of a hypothetical foreground. The dark 
spots mark the locations of the top 5% of foregrounds as measured by the likelihood of the cleaned maps. The lighter, 
smaller spots mark the location of the top 10%, 20% and so forth. The bottom 50% locations are shown as the 
smallest, lightest spots. The top foreground in each case is shown in Table II, along with its location and likelihood. 

It can be seen from Figs. 5-16 that, in most cases, there is a strong clustering of the preferred locations around 
Virgo (and the LSC), its diametrically opposite point, and the region which correspons to Virgo rotated 180° around 
the z-axis. The preferred locations seem to be more scattered for the 3-year Coadded map with both foreground 
models, and for the ILC maps with the HFg-LSC model. 

There is a simple explanation for the spatial distribution seen throughout Figs. 5-16: the HFg's change Prot 
mostly by amplifying and rotating the quadrupole of the original maps. But the quadrupole is even under parity 
transformations n — > —ft, therefore from the point of view of the quadrupole, it is irrelevant if the foreground is at the 
LSC or at its diametrically opposite side. Moreover, since most of the power of the quadrupole lies in its i = 2, m = 
and I = 2, m = ±2 components, a rotation by 180° degrees effects very little change to it. Therefore, if the foreground 
which is distorting the CMB is indeed at the LSC, then its preferred locations in a blind search will degenerate to 
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FIG. 6: Same as before, in the case of the HFg-LSC model applied to the "mask 6" TOH map [TolEll]. 




FIG. 7: Same as before, in the case of the HFg-D model applied to the "mask 6" TOH map [Toll2l|. 

not only the vicinity of the LSC itself, at (I, b) = (284°, 74°), but also to its "dual points" at (I, b) « (104°, -74°) and 
(I, b) w (104°, 74°). This is indeed what seems to happen for most maps - see Table 2. 

It is interesting to notice that, even though the amplitude of the octopole of the HFg-LSC model is subdominant, 
in the HFg-D model there is a substantial change in the octupole after the subtraction. Although the impact of the 
octopole in Prnt is small, there are other statistical tests which are sensitive to it - in particular, the S**- 4 ' 4 ^ statistic 
of Copi et al. 0, which tests for an alignment of the quadrupole and octopole with the ecliptic plane. The 3-year 
Coadded and ILC maps have too high values of S^ 4 ' 4 ) at 99.5% C.L. After subtracting the best-fit Disk foreground 
model, those values come down to 84% C.L. for the Coadded map and to 98% C.L. for the ILC map. 




FIG. 8: Same as before, in the case of the HFg-LSC model applied to the 1-year Coadded map. 
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FIG. 9: Same as before, in the case of the HFg-D model applied to the 1-year Coadded map. 




FIG. 10: Same as before, in the case of the HFg-LSC model applied to the 3-year Coadded map. 

A. Higher multipoles 

If the HFg indeed exists, it may be possible to detect correlations in the higher I components as well. However, 
depending on the model, the spectrum of the foreground can decay with large I, and the increasing number of phases 
means that these correlations will probably difficult to detect. If the foreground is homogeneous (as in the Disk 
model), the higher multipoles can also become important and their presence may affect the alignments between, say, 
I = 2 - 5. 

It is trivial to generalize the tests of Eqs. @-@ for higher multipoles - see, e.g., |2(J |^|. If a given foreground 
cures the quadrupole level and the quadrupole-octopole alignment at the cost of introducing anomalous alignments 
between other multipoles, then the overall likelihood of the resulting map should fall. On the other hand, if the 




FIG. 11: Same as before, in the case of the HFg-D model applied to the 3-year Coadded map. 
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FIG. 13: Same as before, in the case of the HFg-D model applied to the 1-year ILC map. 

original CMB map has another anomalous alignment which is relaxed because of the hypothetical foreground, then 
the likelihood of the resulting map increases even further. 

We have tested the maps of Table II for signs of anomalous alignments between the higher multipoles using an 
unbiased probability with 38 tests, P38 - which includes the level of the quadrupole and the quadrupole-octopole align- 
ment |2 1| . The results are mixed, and probably reveal intrinsic differences between the maps, masks and foreground 
models. For example, whereas the original TOH map with the Mask 6 [Toll2l| has P38 = 3.5 x 10~ 15 , subtraction of 
the best-fit HFg-LSC model leads to an improvement to P38 = 3, 5 x 10~ 10 , while subtraction of the best-fit HFg-D 
model improves it to P38 = 2,4 x 10~ 7 . For the Coadded map (3-year data) with the KP2 mask [l|,|2j; the original 
map has P38 = 4.4 x 10~ 13 , and after subtracting the best-fit HFg-LSC model it only improves to P38 = 1, 1 x 10 -9 , 
while subtraction of the best-fit HFg-D model leads to P38 = 3, 7 x 10~ 8 . Notice that in all cases above, the level of 




FIG. 14: Same as before, in the case of the HFg-LSC model applied to the 3-year ILC map. 
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FIG. 15: Same as before, in the case of the HFg-D model applied to the 3-year ILC map. 



CMB Map - Foreground 


Ptot 


(l,b) 


TOHMask " 


LSC 


0.380 


(136,-51) 


TOHMask " 


Disk 


0.231 


(152,-69) 


TOHMask 6 " 


LSC 


0.488 


(156,-63) 


TOHMask 6 " 


Disk 


0.363 


(159,-73) 


Coadded (lyF 


- LSC 


0.182 


(209,78) 


Coadded (lyr) 


- Disk 


0.144 


(137,50) 


Coadded (3yr' 


- LSC 


0.155 


(314,-10) 


Coadded (3yr) 


- Disk 


0.108 


(122,50) 


ILC (lyr) - 


LSC 


0.219 


(141,77) 


ILC (lyr) - 


Disk 


0.222 


(144,65) 


ILC (3yr) - 


LSC 


0.216 


(133,74) 


ILC (3yr) - 


Disk 


0.155 


(129,61) 



TABLE II: Rotated foregrounds with the highest likelihoods, obtained by maximizing Prot in the cleaned maps. We have 
normalized the foregrounds by their quadrupoles, which we set to be Cf grd ~ 60^K 2 . In terms of the HFg-LSC model, the 
parameter a = 1. In terms of the HFg-D model, the temperature of the 50° x 50° cold spot is Toisk ~ — 130/iK. 

the quadrupole and the quadrupole-octopole alignment alone are responsible for a factor of 10 2 improvement in 
for the TOH map, and for a factor of ~ 10 improvement in P38 for the ILC map. 

V. CONCLUSIONS 

We have presented circumstantial evidence that an extended foreground near the dipole axis could be distorting the 
CMB. The subtraction of such a foreground increases the quadrupole, removes the (anomalous) quadrupole-octopole 
alignment, and dramatically increases the overall likelihood of the CMB maps. Possible physical mechanisms that 
could account for this foreground are the Sunyaev-Zeldovich effect [25j and the Rees-Sciama effect [27|], although it 
should be noted that both options only work in extreme situations that are probably unrealistic. Another possibility 
is that a combination of effects is responsible for the foreground. However, if the Sunyaev-Zeldovich effect due to the 
LSC's gas is indeed responsible for the foreground, it could be directly observed by the Planck satellite |5^| within 
the next few years - see also [5T|. 

We have also shown that the phases of the CMB maps are such that the optimal places for such foregrounds to 
exist would be around the Local Supercluster, its specular image, or the site of the Local Supercluster rotated 180° 
around the galactic polar axis. Furthermore, of the two foreground models analyzed here, the non-uniform foreground 
(HFg-LSC model) seems preferred by the data as it needs a lower overall temperature distortion in order to improve 
the likelihood of the CMB maps. 
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